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Abstract. Interstellar neutral hydrogen atoms flowing into 
the solar system arc attracted by the solar gravitational 
force, repelled by solar hydrogen Ly-a radiation pressure, 
and are ionized, primarily, through charge exchange with 
the solar wind protons. The solar cycle variation of the 
radiation pressure causes the net central solar force to 
fluctuate between attraction and repulsion resulting in the 
modulation of the neutral hydrogen density about the 
usual time independent model. The calculation presented 
here shows that the time dependent downstream density is 
strongly modulated by a large number of travelling neutral 
density waves. The waves possess a continuous range of 
wavelengths as is to be expected for a Maxwellian gas 
subjected to several eleven year solar cycle variations 
during its journey through the solar system. The ampli- 
tudes of the density modulation were found to be quite 
large The backscaltered glow was found to depend on the 
position ol the detector and the phase ol the solar cycle. At 
the most favorable condition a deviation ol the order ol 
25% from the lime independent glow might be observed. 
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1. Introduction 

Neutral hydrogen atoms entering the solar system experi- 
ence a time varying central force resulting from the varia- 
tion of the solar Ly-a radiation pressure. Since the solar 
gravitational force and the radiation force are nearly equal, 
but opposite in direction, the solar cycle variation of 
radiation pressure will make the net central force attractive 
or repulsive depending on the phase of the solar cycle, f he 
neutral atom density along the downstream axis is parti- 
cularly sensitive to the solar cycle cllect. Y or example, in 
the zero temperature limit only those hydrogen atoms 
which experience a net attractive central force while pass- 
ing close to the Sun will reach the downstream axis Thus, 
depending on the hydrogen temperature, the downstream 
time aver need density derived from models not taking into 
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account solar activity variations can differ considerably 
from the instantaneous density. 

A straightforward analytic approach to this time de- 
pendent problem would require that the time dependent 
Boltzmann equation be solved. Fahr ct al. ( 1987) derived a 
formal solution to the Boltzmann equation. It was solved, 
however, only lor the time independent central force 
problem, which is valid for helium but not for hydrogen. 
Recently, Fahr & Scherer (1990) have derived an explicit 
expression for the Courier transformed neutral gas density 
for a time dependent central force. However, this density 
was not back -transformed into real space and time. We 
have taken an alternate and less sophisticated approach by 
following neutral hydrogen atoms through space over 
many solar cvcles for a period of ^ 50 yr. The density is 
obtained by simply counting the number of atoms crossing 
the downstream axis al given time intervals and assigning 
appropriate weights to them. This, however, is only the first 
step. Since we do not perform a full scale simulation and 
have restricted ourselves to following particles with dis- 
crete sots of velocities in the v and y directions to generate a 
provisional density, we have used a houricr transform 
technique (described later) to compute a complete density 
distribution. Our goal in this paper is to study the nominal 
size of the elTect for a time varying radiation force without 
eoinu ihrouuh an expensive and time consuming full scale 
simulation. 


2. Numerical calculation 

Figure 1 shows the geometry of the simulation. The v-axis 
is taken parallel to the neutral hydrogen flow at large 
distances from the Sun. The simulation is simplified by 
assuminu azimuthal symmetry about the .v-axis and is 
carried out in the v r plane. The neutral hydrogen atoms 
are launched Irom a preselected launch space. I he launch 
space is bounded by 220 s. x 5= 20 A Li and 0< i 

It may be thought that launching particles from v - 20 AU 
is too close to the Sun for the Maxwellian distribution to 
hold, however, our result ( Fig. 4) suggests that this approxi- 
mation is adequate. As a further test of our launch space we 
obtain the correct density in the 7=0 K limit, as shown 
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Liter. ) mix is determined empirically, i.e. y ^ increased till 
no particle reaches the downstream axis within 2( )() AU of 
the Sun. All the particles are launched at the same time, in 
1940. The neutral atoms follow trajectories then are non- 
Keplerian because the central force acting on them is time 
dependent. The total radial acceleration of an interplanet- 
ary hydrogen atom is given by 


where is the ratio of the Lvman-alpha radiation pressure 
to the solar gravitational attraction, G is the gravitational 
constant and \1 is the solar mass. The distance of the 
hydrogen atom from the Sun is r and its azimuthal angle 
relative to the upwind axis is <j>. 

The solar cycle variation of the Lyman-alpha radiation 
is described by the time dependence of //. // is known to vary 
from about 0.75 at solar minimum to 1.35 at solar max- 
imum. Based on observations of the solar Lyman-alpha 
radiation we have approximated the variation of // by 


/< 


1 .06 + 0.29 cos 
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where time is given in years AD. 

The density of downwind neutral hydrogen atoms near 
the axis of symmetry is a function of distance to the sun r 
and time /, i.e. den — den(r, M. In order to obtain den(r.r) wc 
have divided the launch area ( Tig 1) into grid points 
spaced by 1/4 AU in the a and y directions. On this space 
grid wc have superimposed a velocity grid, three values of 
V x and V y covering the range appropriate to a Maxwellian 
velocity distribution ot an 8000 K gas. From each point of 
the combined space velocity grid one particle is launched. 
The orbits of all particles are lollowed and those which 
cross the downwind region are put into various bins 
according to the position r and the time t of their down- 
wind axis crossing. Fach bin covers a spatial region of 
dimension A.v Ay centered about the position r = x and 
y — 0. The particles are collected over a time interval A / at 
time r. A.v and Ar were both taken to be I AU and Ai was 
taken to be one year. After all the atoms passed through the 
downwind region, the density denfr, t) was found by adding 
ail properly weighted particles w hich passed through each 
bin. The weighting factor lakes into account all of the 
following points: 

(a) The probability that the atom at the launch posi- 
tion has the velocity components I A and V y 

(b) The volume of the launch space, which is propor- 
tional to 2nv lfl , where y in is the y-coordinatc at the launch 
point. This factor arises from our assumption of azimuthal 
symmetry. 

(c) The surviv al probability of the atom against ioniza- 
tion. 

|d) The time span during w hich a particle remains in its 
particular cell in t he downwind region. 


We will briefly discuss the various terms of the weighting 
factor. 

We assume the neutrals in our launch space to have a 
temperature / and therefore a Maxwellian velocity dis- 
tribution centered on V x — L h and l { 0, where I h is the 

bulk velocity of the neutrals taken lo be 20kms The 
density of the neutrals in the launch region is taken to be 
uniform and the velocity distribution of the neutrals in the 
launch region is proportional to 


exp 


2A7' <,/ -^ +K? 


where m is the hydrogen atom mass, / the temperature and 
k the Boltzmann constant. 

The downwind neutral density is proportional to the 
survival probability surv(/\ /) of neutrals against inter- 
planetary ionization processes. The predominant ioniza- 
tion process is charge exchange with solar wind protons, 
and the probability of survival against this process is given 
by exp( fd//r) = cxp( — J« p fT F rc , dr). Here r is the lifetime 
of a H atom against ionization by charge exchange. /? p is 
the density of protons at (/*, 0), a is the charge exchange 
cross section, and F rc , is the relative velocity between the 
neutral atoms and protons. The angular momentum of a 
particle at any point along a trajectory is constant and 
given by / = nir 2 dr/i/dr Hence, we can substitute dr 
= {mr 2 l) d(j> into the expression for surv. The final ex- 
pression for surv is exp( - f F rc ,/nr d</> /). 

The time interval during which a particular particle 
remains in a given bin on the downstream axis is propor- 
tional to Ay/r v . Particles which traverse more than one bin 
along the axis are in this approximation placed in the first 
bin which they enter. 

Based on these considerations we obtain dcn(M) by 
summing over all the weighted particles in the bin defined 
by r and I: 
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3. Results and discussion 

The zero order density calculated b\ this method uses a 
limited set of velocities. So wc use Lourier transform 
method to compute the complete density as discussed 
below. We computed the zero order neutral hydrogen 
atom density for the downstream region lor various years 
between 1972 to 1992 by counting the particles collected in 
the various bins. The statistical accuracy of the provisional 
density varied with the distance from die Sun. The accur- 
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acy was of high order <5% at heliocentric distance 
> 100 AU. Near the Sun ( < 50 AIJ) the accuracy declines 
to 10% or more. Since it is clear that the neutral hydrogen 
atom density will carry information about the solar cycle 
variation of the radiation pressure we have calculated the 
power spectrum of the density fluctuations as a function of 
the wave vector k and frequency < The Fourier integral 
transform of den(r. i) is given by 


G(fc,<u) 



den0\ /) e ,ar vm dr dr. 


(4) 


The power spectrum in ru. }(i(k.to)\ 2 , for a particular value 
of k is shown in Fig. 2. The first peak at is the time 

independent part of the density while the peak at 
(l) ^2n 11 yr 1 corresponds to the II yr period arising 


Fig. 1. The launch geometry is shown. The shaded 
area is the region from which the atoms are laun- 
ched. The black dots schematically represent the 
discrete launch points 


Fig. 2. The power spectrum of the full 
< density |( i[k % rv)| 2 for A 0.1 Al) 1 is plot- 
ted against <^(yr '), 1 his calculation was 
■ performed for a 20 yr period from 1972 to 
1992 

1 1.2 


naturally from the solar activity cycle. Since we did not 
carry out an exact simulation it is not surprising that the 
expected peak is not exactly at m-0. I he third peak 
may be an artifact of the approximation method used here. 
This result has enabled us to write the density dcn(r, f) as 
denlr. U dcn tl (r)[ 1 + </,(r)cxp( - where dcn 0 (r) is the 
time independent density. The absence of neutral hydrogen 
atoms close to the Sun (r-0) imposes the condition that 
den o (0} 0. T he time independent density den 0 (r) is of the 
form u(l e hr ). The time independent density is of the 
exponential form if the solar gravitational force and radi- 
ation force cancel each other (Axford 1972; Iscnberg 1986). 
The constants a and h w'ere obtained by a least squares fit 
to the simulated time dependent data. The fluctuating part 
r/j(r)e\pi umM was obtained from the ratio [den(M) 
den,, { r) ] den„(r). 
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The Fourier transform of the oscillating density is given 
by 


i}(k % v>) = 


iln) 2 



d j(r)e ’* Jnl e tikr ""drd. 


(2n) 

1 

(2tt) 


dr 


d/ 


d x {r)c ' kr drdia ) 0 - ro) </(A’) m). 


( 5 ) 


The power spectrum in k space for the oscillating density, 
|<y(A)| 2 , computed lor /=1988 is shown in Pig. 3. 1 he 
expression for (j(k) was obtained from I q. (5). Lor the year 
1988 there are three clear peaks in k space. These three 
peaks correspond approximately to the three (FJ velocit- 
ies of launch. An exact simulation using a finer velocity grid 
would reveal a continuous spectrum in k. This is a reflec- 
tion of the Maxwellian distribution of velocities. There is 
also an interesting feature in the power spectrum in a). The 
density power spectrum is in principle a ^-function in 
(u-spacc because all the atoms were subjected to the same 
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i Lig. 3. The power spectrum of the re- 

| sidual neutral density \(}[k)\ 2 is plotted 

; against wave vector k { AU l ), where </|A) 
= (e lU d,(r) dr. This calculation was per- 
formed over a distance between 5 to 
■-* 140 At 1 for the year 19S8 
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j Fig. 4. The last square fitted neutral den- 
sities for 1988 (dots) and the modulated 
| density (solid line) are plotted against the 
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I I yr variation ol the net central force. I he finite spatial 
and temporal simulation of the neutral atom flow through 
the solar system, however, “broadens" the ^-function 
which shows up as the broad peak in Fig. 2 . 

Since there is a continuous spectrum in k the density 
should be written in a Fourier integral from: 


dcn(rj) = dcn 0 (r) 




0(Jc,w)c Hfcr 


n dkd(o 


= dcn 0 (r) 



i/(k) 
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= den n (r) 


1 +c 


c)(k)c ikr 



( 6 ) 


We have accordingly used Ecp (6) to calculate the total 
density in 1988 and the results are shown in Fig. 4 . Of 
course, the physical density is the real part of dcn(r, r). It is 
interesting to note that the size of the modulation is quite 
large in the downstream direction. This is consistent with 
the limiting case of a cold gas, T= 0 K (Fig. 5 ), where the 
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Fig. 5 . The downstream neutral hydrogen 
density for a 0 K temperature gas as ob- 
tained from the simulation is plotted 
against the heliocentric distance 
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Fig. 6. The quantity \{n r 1 ) dr, which is 
proportional to the backscattered inten- 
sity, is plotted for the year 1988 against 
heliocentric distance for cases (i) the lime 
independent density (dots) and (ii) modul- 
ated density (solid line) 
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modulation is 1 00° n on the downwind axis giving rise to 
travelling hydrogen rich and hydrogen empty regions 

Would it be possible to see these travelling waves by 
using the IJ V photodetectors aboard the deep space probes 
Pioneers 10/1 1 and Voyagers 1/2? We have carried out a 
demonstrative calculation to explore the cfTccl of the solar 
cycle driven wave on the backscattered intensity We 
calculated J(n/r 2 )dr (a quantity proportional to intensity in 
the single scattering approximation) where n is the hydro- 
gen density and r is the heliocentric distance to the 
scattering volume. For simplicity we assumed an isotropic 
phase function and a radial field of view. We calculated the 
backscattered hydrogen Ly-ot intensity for both time de- 
pendent and time independent density for the year 1988, i.e. 

(i) n = dcn 0 (r) and 

(ii) n — dcn 0 (r, /)[! +di(r)exp( — iatf)]. (7) 

In Fig. 6 wc have plotted the backscattered intensities in 
the downwind direction for cases (i) and (ii). Any deviation 
in the backscattered glow due to the travelling wave is seen 
to depend on the spacecraft position and the time of 
observation and can be as large as ^25%. T he Pioneer 10 
glow data may contain evidence of this effect since it is 
moving in the downstream direction where the density 
variations arc significant and in principle observable. The 
density modulations might also be large enough to be 
directly detected by a neutral particle in-situ detector 


4. Summary 

The effect of solar cycle radiation pressure variations on 
the downstream neutral hydrogen density has been ex- 
plored. The Fourier transform of the density in k and a> 
space was found to depend on a spectrum of k values but 
only on a single o)( — 27t/l 1 yr s 1 ). There is clear evidence 
for a series of travelling waves with a continuous spread of 
wave vectors. Thus, a spacecraft going downstream should 
be able to detect a deviation in the backscattered glow 
intensity, relative to that expected on the basis of a time 
independent density, for selected detector positions and 
observation times. 
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